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Abstract 

We investigate qualitatively and quantitatively the impact of the 
general relativistic gravito-electromagnetic forces on hyperbolic orbits 
around a massive spinning body. The gravito-magnetic field, which 
is the cause of the well known Lense-Thirring precessions of elliptic 
orbits, is generated by the spin S of the central body. It deflects and 
displaces the trajectories differently according to the mutual orienta- 
tion of S and the orbital angular momentum L of the test particle. The 
gravito-electric force, which induces the Einstein precession of the per- 
ihelion of the orbit of Mercury, always deflects the trajectories inward 
irrespective of the L S orientation. We numerically compute their 
effect on the range r, radial and transverse components v r and v T of the 
velocity and speed v of the NEAR spacecraft at its closest approach 
with the Earth in January 1998 when it experienced an anomalous 
increase of its asymptotic outgoing velocity Voo o of 13.46 ± 0.01 mm 
sec -1 ; while the gravito-electric force was modeled in the software used 
to process the NEAR data, this was not done for the gravito-magnetic 
one. The range-rate and the speed are affected by general relativistic 
gravito-electromagnetism at 10~ 2 (gravito — electric) — 10~ 5 (gravito- 
magnetic) mm sec -1 level. The changes in the range are of the order 
of 10 -2 (gravito — magnetic) — 10 1 (gravito — electric) mm. 



PACS: 04.80.Cc; 95.10.Ce; 95.55.Pe 



1 Introduction 

In this paper we will investigate the effects of general relativity, in its weak- 
field and slow-motion approximation, on unbound, hyperbolic orbits of test 
particles approaching a body of mass M and angular momentum S. We 
will consider both the gravito-magnetic and the gravito-electric relativistic 
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Figure 1: Essential features of the unperturbed Newtonian hyperbola. P is 
the position of a test particle at time t with respect to the body located at the 
origin of the chosen reference frame. The smallest distance is r m [ n = o(e— 1). 
The angle between r and v is <f>. M here is not to be confused with the central 
body of mass M. 

forces induced by the stationary and static components, respectively, of the 
field of M [1]. 

The results obtained will be applied to some realistic planet-spacecraft 
scenarios in the Solar System to see if the predicted effects fall within the 
current or future sensitivity level of the ranging techniques. In particular, we 
will examine the so-called flyby anomaly consisting of unexplained changes 
of the asymptotic outgoing velocities of some spacecraft (Galileo, NEAR, 
Cassini and MESSENGER) that occurred at their closest approaches with 
the Earth [2, 3]. 

Let us recall some basics of the Newtonian hyperbolic orbit [4] which 
represents, in this case, the reference unperturbed path; see Figure 1. Its 
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equation is 

1 + e cos / 

where p = a(e 2 — 1) is the semilatus rectum, and / is the true anomaly 
considered positive in the anti-clockwise direction from the point of closest 
approach A; a is the semi-major axis and e > 1. The speed is 



(2) 




so that 



^. (3) 

a 

The asymptotic ingoing and outgoing speeds, which are equal to Voo in the 
unperturbed case, are denoted with ; and Q , respectively. The angle 
<j> between the r and v is defined by 

. 1 + ecos/ 

sinc£>= — . , (4) 

y/l + e 2 + 2ecos/ 

& sin j~ 

cosd> = (5) 

Vl + e 2 + 2ecos/ 

For x = 0, y = ±p, i.e. f = f = ±vr/2, it is 

v x0 = vs'mcf) = (I— , (6) 




v y o = vcos(j) = Te\ • (7) 



2 The gravito-magnetic force 

In this section we will deal with the effect of the general relativistic gravito- 
magnetic force on the hyperbolic motion of a test particle approaching a 
spinning body of mass M and angular momentum S. 

Let us briefly recall that, in the weak-field and slow-motion linear ap- 
proximation of general relativity, the off-diagonal components goi,i = 1)2,3 
of the space-time metric tensor, related to the mass-energy currents of the 
source, induce a gravito-magnetic field B g [1] by analogy with the magnetic 
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field caused by moving electric charges in the linear Maxwellian electro- 
magnetism. Far from an isolated rotating body, the gravito-magnetic field 
becomes [5] 

B g = ~[5-3(5-f)f], (8) 

where G is the Newtonian gravitational constant and c is the speed of light 
in vacuum. It exerts the non-central Lorentz-like acceleration [1] 

A GM = --v X B g (9) 

upon a test particle moving with velocity v. For ordinary astronomical 
bodies like, e.g., the Earth and the Sun, ^4 GM is many orders of magnitude 
smaller than the Newtonian monopole A N = GM/r 2 , so that it can be 
considered small perturbation. 

The action of eq. (9) in the case of unperturbed close orbits, giving 
rise to, e.g., the Lense-Thirring precession of the ellipse of a test particle 
[6, 7, 8], has been the subject of intense activity, both theoretically and 
observationally, in recent times [9]. 

Here we will consider as reference path of a Newtonian hyperbolic tra- 
jectory. In order to work out the effects of the gravito-magnetic field on it, 
we will numerically integrate the equations of motion in cartesian rectangu- 
lar coordinates [10] for some particular orbital geometries over a time span 
including the epoch of closest approach to M which is assumed located at 
the origin of the coordinate system. 

2.1 Qualitative features for equatorial and polar osculating 
orbits 

First, we will consider a trajectory lying in the equatorial plane of the ro- 
tating body for the cases of co-rotation (Figure 2) and counter-rotation 
(Figure 3) of the particle's radius vector r with respect to the diurnal ro- 
tation of M whose spin S is assumed to be directed along the positive z 
axis (anti-clockwise diurnal rotation). Just for illustrative purposes, we as- 
sumed the Earth as source of the gravitational field and re-scaled its gravito- 
magnetic force by a factor 10 10 in such a way to still keeping the condition 
A GM /A^ < 1 valid over the entire orbit. It turns out that the perturbed 
trajectory remains confined in the equatorial plane of the central body; for 
co-rotation the path is deflected outward with respect to the unpertubed 
hyperbola, while for counter-rotation it is deflected inward. Indeed, for 
equatorial orbits A GM is entirely in-plane because B g is directed along the 



4 




Figure 2: Effect of the gravito- magnetic force on the hyperbolic motion of 
a test particle around an astronomical rotating body located at the origin 
of the depicted frame. The body's spin S is directed along the positive z 
axis, i.e outside the figure. Red line: unperturbed hyperbola. Blue dashed 
line: perturbed orbit. For illustrative purposes we choose the Earth as 
central body and re-scaled the magnitude of its gravito-magnetic force by 
10 10 so that A GM -/A N = 0.4 at perigee. We adopted the initial conditions 
xo = 0, y = -p = -a(e 2 - 1), z = 0, v x0 > 0, v y0 > 0, v z0 = 0; the 
particle moves in the equatorial plane of the spinning Earth in such a way 
that the radius vector rotates in the same sense with respect to the Earth, 
i.e anticlockwise. We used a = 8493.326 km, e = 1.81. The perturbed orbit 
is deflected outward with respect to the unperturbed one. 
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Figure 3: Effect of the gravito-magnetic force on the hyperbolic motion of 
a test particle around an astronomical rotating body located at the origin 
of the depicted frame. The body's spin S is directed along the positive z 
axis, i.e outside the figure. Red line: unperturbed hyperbola. Blue dashed 
line: perturbed orbit. For illustrative purposes we choose the Earth as 
central body and re-scaled the magnitude of its gravito-magnetic force by 
10 10 so that A GM /A^ = 0.4 at perigee. We adopted the initial conditions 
xq = 0, yo = p = a(e 2 — 1), z$ = 0, v x $ > 0, v y o < 0, v z q = 0; the particle 
moves in the equatorial plane of the spinning Earth in such a way that the 
radius vector rotates in the opposite sense with respect to the Earth, i.e. 
clockwise. We used a = 8493.326 km, e = 1.81. The perturbed orbit is 
deflected inward with respect to the Unperturbed one. 



negative z axis; for co-rotating particles it is radially directed outward at 
the point of closest approach and decreases the gravitational pull felt by the 
orbiter, while for counter-rotating probes it is radially directed inward at the 
pericentre and increases the net gravitational acceleration. The flyby epoch 
is left almost unaffected. By taking the difference between the integrated 
perturbed and unperturbed orbits sharing the same initial conditions, it 
can be shown that in the co-rotating case the radial velocity, which is one 
of the direct observables in real planet-spacecraft close encounters, experi- 
ences an increase with respect to the unperturbed one just around the flyby 
epoch, while the radial components of v<x> \ and Vqq are left unaffected. In 
the counter-rotating case v r decreases at the closest approach. Concerning 
the body-centric range r, it turns out that it suddenly increases (decreases) 
around the flyby epoch for the co-(counter-)rotating case and remains about 
at that level also after the flyby. 

Let us, now, consider the case in which the unperturbed hyperbola en- 
tirely lies in an azimuthal plane, e.g. the {yz} plane, containing the spin 
S of the central body. Now, since v is contained in the same plane of B g , 
the gravito-magnetic acceleration is out-of-plane, so that it can be expected 
that the perturbed trajectory will be displaced along the x axis. This fact is 
confirmed by a numerical integrations shown in Figure 4 and Figure 5 from 
which the displacement of the orbit perpendicularly to the initial osculating 
plane is apparent. 

2.2 The flyby anomaly: the NEAR case 

Many spacecraft launched in interplanetary missions make use of one or 
more Earth's flybys in order to gain or lose the heliocentric energy required 
to reach their far targets (planets, asteroids, comets) without using huge 
and expensive amounts of propellant [11, 12]. In the case of GALILEO 
(twice), NEAR, Cassini and MESSENGER unexplained variations in Voo 
were detected [2, 3]; the largest one was measured at the close encounter of 
NEAR with the Earth that occurred in 1998 and amounts to 



The unperturbed hyperbola of NEAR is depicted in Figure 6. 

Anderson et al. in [3] derived an empirical formula which satisfactorily 
fit all the six flyby anomalies measured so far. It is 



Avon = 13.46 ± 0.01 mm sec 



-l 



(10) 




(cos (5; - cos^ ) = (3.099 x 10~ 6 ) (cos 5; - cos5 ) , (11) 
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Figure 4: Effect of the gravito-magnetic force on the hyperbolic motion of 
a test particle around an astronomical rotating body located at the origin 
of the depicted frame. The body's spin S is directed along the positive z 
axis, i.e outside the figure. For illustrative purposes we choose the Earth as 
central body and re-scaled the magnitude of its gravito-magnetic force by 
10 10 so that A GM /A^ = 0.4 at perigee. We adopted the initial conditions 
= 0, yo = -p = -a(e 2 - 1), z = 0, v x0 = 0, v y0 > 0, v z0 > to have 
the spacecraft initially moving in the osculating {yz} plane. The perturbed 
trajectory is displaced along the x axis on the upper left corner. 
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Figure 5: Effect of the gravito- magnetic force on the hyperbolic motion of a 
test particle around an astronomical rotating body located at the origin of 
the frame shown with its spin directed along the positive z axis, i.e outside 
the figure. For illustrative purposes we choose the Earth as central body 
and re-scaled the magnitude of its gravito-magnetic force by 10 10 so that 
A M /A^ = 0.4 at perigee. We adopted the initial conditions xq = 0, yo = 
p = a(e 2 — 1), zo = 0, v x q = 0, v y o < 0, v z $ > to have the spacecraft 
initially moving in the osculating {yz} plane. The perturbed trajectory is 
displaced along the x axis on the lower right corner. 
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Figure 6: Unperturbed hyperbola of NEAR; its osculating plane is tilted by 
108 deg to the {x, y} plane assumed coincident with the Earth's equator. 
The starting point is in the right upper corner (xo > 0, yo > 0, zq > 0). 
The Earth has to be imagined located at the origin of the coordinate system 
with its spin S directed along the positive z axis. 
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where oj and R are the Earth's angular speed and equatorial radius, respec- 
tively, while 5{ and 5 Q are the ingoing and ongoing geocentric declinations, 
respectively. Concerning possible explanations in terms of known physics, 
many dynamical effects (tides, atmospheric drag, Earth's albedo, solar wind, 
terrestrial magnetic field, spin-rotation coupling between electromagnetic 
waves and spacecraft and Earth rotations) were preliminarily considered by 
Lammerzahl et al. in [13]; an order-of-magnitude approach was followed by 
confronting the magnitude of the accelerations induced by standard forces 
considered with the one which is assumed to be responsible of the flyby 
anomaly, i.e. ~ 10 -4 m sec -2 . As a consequence, all the effects investigated 
were discarded. However, we note that detailed analyses are in order: in- 
deed, even if some dynamical effect, standard or not, was found to produce 
an acceleration with the right order of magnitude, it might happen that its 
signature on the observable quantities is not correct; that is, it may induce, 
for instance, a decrease of the radial velocity. Mbelek in [14] suggested that 
the special relativistic Doppler effect may explain the formula by Anderson 
et al. [3]. Among various explanations in terms of non-conventional physics 
put forth so far, McCulloch in [15] proposed a mechanism based on the 
hypothesis that inertia is due to a form of Unruh radiation and varies with 
acceleration due to a Hubble-scale Casimir effect. It qualitatively reproduces 
the latitude-dependence of eq. (11) and is quantitatively in agreement with 
three of the six measured flybys. 

Since eq. (11) contains a term including quantities like the speed of light 
and the first power of the Earth's angular speed which enters general rela- 
tivistic gravito-magnetic effects, it seems interesting to apply our previous 
results concerning the influence of B g on hyperbolic orbits to the NEAR's 
flyby. Note that the gravito-magnetic force was not modeled in the software 
used to process the NEAR data. By using the state vector of NEAR at the 
flyby epoch (J D Anderson, private communication to the author, November 
2008) referred to a geocentric equatorial frame {a;, y, z}, and eq. (9) it turns 
out that at the point of closest approach to Earth along the flyby trajectory 

A™ = 3.3 x 10" 10 m sec" 2 , (12) 

A™ = 7.5 x 10~ n m sec" 2 , (13) 

= -1.7 x 10" 10 m sec" 2 , (14) 

so that 

i GM = 3.8 x 10~ 10 m sec" 2 . (15) 

We will now use a numerical integration of the equations of motion per- 
turbed by eq. (9). We look at a time span starting from the flyby epoch 
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Figure 7: Change in the radial velocity Av r induced by the gravito-magnetic 
force on NEAR at the Earth's flyby (1998-Jan-23 07:22:56 Coordinate Time 
(CT)). It has been obtained by taking the difference between the integrated 
perturbed and unperturbed trajectories sharing the same initial conditions 
of Table 1. They have been obtained with the HORIZONS software by 
NASA JPL at 1998-Jan-23 07:00:00 CT and correspond to an instant 1353 
sec before the flyby. Reference frame: ICRF/J2000.0. Coordinate system: 
Earth Mean Equator and Equinox of Reference Epoch. The maximum effect 
Av™ ax = —5 x 10~ 5 mm sec -1 occurred just at the flyby. 



and extending in the future after it: the chosen initial conditions are in 
Table 1. We will consider the changes of the velocity along the radial r , 
transverse f and out-of-plane v directions; v is directed along the orbital 
angular momentum, while f = v X f is not directed, in general, along v. 

The results for Av r , Av T , Av and Ar are shown in Figure 7-Figure 
10, respectively. They have been obtained by subtracting the unperturbed 
orbit from the perturbed one, both numerically integrated with the initial 
conditions of Table 1. 

The gravito-magnetic force of the Earth decreased the radial velocity of 
NEAR by ~ 10 -5 mm sec -1 just at the flyby, while the transverse velocity 
was augmented after the flyby up to 10 -5 mm sec -1 level; it turns out that 
the normal velocity was affected at an even smaller level. The total speed 
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Table 1: Initial conditions used for NEAR obtained with the HORIZONS software by NASA, JPL at 1998-Jan-23 
07:00:00 CT (Coordinate Time, defined as the uniform time scale and independent variable of the ephemerides 
http://ssd.jpl.nasa.gOv/7horizons_doc#timesys), i.e. 1353 sec before the flyby. Reference frame: ICRF/J2000.0. 
Coordinate system: Earth Mean Equator and Equinox of Reference Epoch. 

xq (km) i/Q (km) zq (km) v x q (km sec" 1 ) v y o (km sec" 1 ) v z q (km sec -1 ) 



4,496.885594909381 6,930.477153733549 13,199.11503591246 -1.712684317202157 -8.679677119077454 -4.455285829060190 



Av r , mm/ sec 




Figure 8: Change in the transverse velocity Av T induced by the gravito- 
magnetic force on NEAR at the Earth's flyby (1998-Jan-23 07:22:56 CT). 
The maximum value reached is of the order of 10~ 5 mm sec -1 . 
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, sec 



Figure 9: Change in the speed Av induced by the gravito-magnetic force 
on NEAR at the Earth's flyby (1998-Jan-23 07:22:56 CT). The maximum 
effect Au max = 2 x 10~ 5 mm sec -1 occurred just at the flyby while Voo is 
left almost unchanged. 
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Figure 10: Change in the radial distance Ar induced by the gravito-magnetic 
force on NEAR at the Earth's flyby (1998-Jan-23 07:22:56 CT). The maxi- 
mum effect Ar max = —6 x 10 -2 mm took place after the flyby. 
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v was increased up to 10 mm sec - i.e. six orders of magnitude smaller 
than the observed increment. The geocentric range of the spacecraft was 
reduced by about 10~ 2 mm. 



3 The gravito-electric force 

In the Post-Newtonian approximation of order 0(c~ 2 ), the acceleration in- 
duced by the gravito-electric component of the field of a static mass M is, 
in standard isotropic coordinates, [10] 



A 



GE 



GM 



4GM v 2 4(r-v) 

V -v -\ V 

y4 ^->3 ^"3 



(16) 



Concerning the influence of eq. (16) on the hyperbolic motion, it has only 
radial and transverse components, so that no departures from the osculating 
plane occurs. From Figure 11 it turns out that the test particle is deflected 
inward with respect to the unperturbed hyperbola. 

The gravito-electric acceleration experienced by NEAR at the point of 
closest approach to earth along the flyby trajectory was 

A^ E = 9.5 x 10" 10 m sec" 2 , (17) 

A^ E = -5.26 x 10" 9 m sec" 2 , (18) 

A^ E = 3.42 x 10" 9 m sec" 2 , (19) 

so that 

,4 GE = 6.35 x 10" 9 m sec" 2 . (20) 

The impact of the gravito-electric force 1 on v r , v T , v and r of NEAR 
at its flyby are depicted in Figure 12-Figure 15, respectively; their patterns 
are quite different from the gravito-magnetic ones; the outgoing asymptotic 
velocity is changed by an amount of the order of 10~ 3 mm sec" 1 . The 
maximum variations of the range rate, the transverse velocity, the speed 
and the range are of the order of 10~ 2 mm sec" 1 and 10 1 mm, respectively. 



1 Contrary to the gravito-magnetic one, it was modeled in the software used for pro- 
cessing the data. 
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Figure 11: Effect of the gravito-electric force on the hyperbolic motion of a 
test particle around an astronomical body located at the origin of the frame 
shown. Red line: unperturbed hyperbola. Blue dashed line: perturbed orbit. 
For illustrative purposes we choose the Earth as central body and re-scaled 
the magnitude of its gravito-electric force by 10 8 so that A GE /A^ = 0.1 at 
perigee. We adopted the initial conditions xq = 0, yo = —p = —a(e 2 — 
1), z = 0, v x0 > 0, v y0 > 0, VgQ = 0. We used a = 8493.326 km, e = 1.81. 
The perturbed orbit is deflected inward with respect to the unperturbed 
one. 
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Figure 12: Change in the radial velocity Av r induced by the gravito-electric 
force on NEAR at the Earth's flyby (1998-Jan-23 07:22:56 CT). It is the 
difference between the integrated perturbed and unperturbed trajectories 
sharing the same initial conditions of Table 1. They have been obtained 
with the HORIZONS software by NASA JPL at 1998-Jan-23 07:00:00 CT 
and correspond to an instant 1353 sec before the flyby Reference frame: 
ICRF/J2000.0. Coordinate system: Earth Mean Equator and Equinox of 
Reference Epoch. 
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Figure 13: Change in the transverse velocity Av T induced by the gravito- 
electric force on NEAR at the Earth's flyby (1998-Jan-23 07:22:56 CT). 



Av, mm/sec 




, sec 



Figure 14: Change in the speed Av induced by the gravito-electric force on 
NEAR at the Earth's flyby (1998-Jan-23 07:22:56 CT). 
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Figure 15: Change in the radial distance Ar induced by the gravito-electric 
force on NEAR at the Earth's fiyby (1998-Jan-23 07:22:56 CT). 

4 Discussions and conclusions 

We investigated qualitatively and quantitatively the impact of general rel- 
ativity, in its weak-field and slow-motion approximation, on unbound hy- 
perbolic orbits around a massive spinning body. We considered both the 
gravito-magnetic and the gravito-electric terms; the first one is responsi- 
ble of the Lense-Thirring precessions of elliptic orbits, while the second one 
causes the well known Einstein precession of the perihelion of Mercury of 
43.98 arcsec cty -1 . The gravito-magnetic force deflects an equatorial trajec- 
tory inward or outward with respect to the unperturbed hyperbola according 
to the mutual orientation of the orbital angular momentum L of the particle 
with respect to the spin S of the central body. For osculating orbits lying 
in a plane which contains S there is also a displacement in the out-of-plane 
direction. The gravito-electric force is not sensitive to the L — S orien- 
tation and deflects the trajectory inward with respect to the unperturbed 
hyperbola. 

We applied our results to the flyby anomaly experienced by the NEAR 
spacecraft at its close encounter with the Earth on January 1998 when its 
asymptotic outgoing velocity Voo was found larger than the ingoing one by 
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13.46 ± 0.01 mm sec ; contrary to the gravito-electric force, the gravito- 
magnetic one was not modeled in the software used to process the NEAR 
data. From numerical integrations of the perturbed equations of motion in 
a geocentric equatorial frame with rectangular cartesian coordinates over a 
time span extending in the future after the flyby epoch, we quantitatively 
investigated the changes in the radial and transverse components of the 
velocity v r and v T , the speed v and the range r of NEAR induced by the 
general relativistic gravito-electromagnetic forces. Concerning the range, its 
variations are at the 10 -2 level for the gravito-magnetic force and 10 1 mm 
level for the gravito-electric one. The radial and transverse velocities and 
the speed are affected at 10 -2 (gravito-electric) — 10 -5 (gravito-magnetic) 
mm sec -1 level. 
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